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We study the thermodynamics of large N pure 2+1 dimensional Yang-Mills theory on a 
small spatial S 2 . By studying the effective action for the Polyakov loop order parameter, 
we show analytically that the theory has a second order deconfinement transition to a phase 
where the eigenvalue distribution of the Polyakov loop is non-uniform but still spread over 
the whole unit circle. At a higher temperature, the eigenvalue distribution develops a 
gap, via an additional third-order phase transition. We discuss possible forms of the full 
phase diagram as a function of temperature and sphere radius. Our results together with 
extrapolation of lattice results relevant to the large volume limit imply the existence of a 
critical radius in the phase diagram at which the deconfinement transition switches from 
second order to first order. We show that the point at the critical radius and temperature 
can be either a tricritical point with universal behavior or a triple point separating three 
distinct phases. 
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1. Introduction 

In this note, we follow to study the thermodynamics of large N pure 2+1 dimen- 
sional Yang-Mills theory on a spatial S 2 with radius much smaller than the scale set by 
the dimensionful coupling of the gauge theory. In this limit, the dimensionless coupling 
A = g 2 NR is small, so the thermodynamics can be studied in perturbation theory. The 
thermal partition function is given by the path integral for the Euclidean theory on S 2 x S , 
where the S 1 has radius given by the inverse temperature (3 = T _1 . From this path integral 
expression, we integrate out all modes apart from the trace of the Polyakov loop,0 

u = ^tr(U) = ^trPe l ^ A , 

giving us an effective action for u, the standard order parameter for confinement. This 
effective action takes the form (valid for \u\ < 1/2) 

S eff (u) = f(T, X)\u\ 2 + X 2 b(T)\u\ 4 + 0(X 4 ) . (1.1) 

For low temperatures, / is positive, and so the saddle-point configuration has u = 0. At 
some temperature T# = + 0(A) , / becomes negative, so u = is no longer the 
minimum action configuration. The three-loop calculation in this paper shows that the 
coefficient b is positive at the critical temperature, so as / becomes negative, \u\ develops 
an expectation value gradually, resulting in a second-order phase transition. 

The situation here is in contrast to the 3+1 dimensional case, where the analogous 
calculation [jl[ revealed a negative value for b{Tn)- In that case, the deconfinement tran- 
sition is first order, characterized by a discontinuous change in the eigenvalue distribution 
for the unitary matrix U from the uniform distribution (eigenvalues equally distributed 
around the unit circle) to a non-uniform distribution with a gap (i.e. a region of the circle 
where no eigenvalues are present). In terms of the eigenvalue distribution, the continuous 
transition we find here corresponds to a continuous change from the uniform distribution 
to an increasingly non-uniform distribution. As argued in || , at a temperature 

T 2 = T H + ^X 2 b(T H )/f'(T H ) 

the minimum value of the eigenvalue density will reach zero, as depicted in figure 1, and 
the eigenvalue distribution will develop a gap for higher temperatures. This results in 

1 Here A is averaged over the S 2 . 
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Fig. 1: Distribution of Polyakov loop eigenvalues on unit circle (horizontal 
axis) in confined phase (a), in gapless phase above second order deconfinement 
transition (b), at third-order gapping transition (c), and in high-temperature 
gapped phase (d). 

an additional third-order phase transition of Gross- Witten type, so the phase diagram 
for pure Yang-Mills theory on S 2 x S 1 contains at least three distinct phases (uniform, 
non-uniform, gapped). 

The basic setup for our calculation and the calculation itself are presented in sec- 
tions 2 and 3 of this paper. We then consider the implications of our results for the full 
phase diagram as a function of temperature and spatial radius. In section 4, we argue 
that in the high-temperature limit for fixed radius the theory effectively reduces to pure 
two-dimensional Yang-Mills theory on a spatial S 2 . This theory has a third order phase 
transition as the radius of the sphere is varied 0, suggesting that our phase diagram has 
an additional third order transition line coming from large temperature along the curve 

TR ~ 1/(XR) . 

It is tempting to conjecture that this third-order transition line connects up with the one 
emerging from the critical temperature at small radius (as in figure 6b). Indeed, these are 
qualitatively very similar: while the latter transition is associated with the development of 
a gap in the eigenvalue distribution for the Polyakov loop, the two-dimensional Yang-Mills 
theory transition is associated with the development of a gap in the eigenvalue distribution 
for the Wilson loop around the equator of the sphere (or any other maximal-area non- 
intersecting loop) 0. On the other hand, we have not been able to show any direct relation 
between the two order parameters (which we can choose to be the minimum value of 
the eigenvalue density for these two Wilson loops). Further, it is possible that the high 
temperature transition ceases to be sharp as we decrease the temperature. 

In section 5, we consider the fate of the actual deconfinement transition line as the 
radius is increased and the possible forms for the full phase diagram. At large volume, 
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we can appeal to lattice results, which for finite iV suggest that there is a second order 
deconfinement transition for iV = 2, 3 (and possibly N = 4), a weak first order transition 
for N = 5 and a stronger first-order transition for iV = 6 0. It is believed that the 
transition should continue to strengthen as iV is increased. If this is correct, the large 
iV theory would be expected to exhibit a first order deconfinement transition, and there 
must be a point along the deconfinement transition line at some critical radius, where the 
deconfinement transition switches from second order to first order. We argue that there are 
two qualitatively different behaviors of the Polyakov loop effective potential that can lead 
to such a point. These correspond either to a tricritical point at which both quadratic and 
quartic terms in an effective action vanish, or to a triple point, separating three distinct 
phases. The simplest possible phase diagrams consistent with the available information 
are presented in figure 6. 

We conclude in section 6 with a few comments on the implications for a possible 
gravitational dual theory of the existence of three distinct phases. 

2. The set up 

In this section we will briefly review the results and techniques developed in JI| , |J , || 
and apply them to analyze the thermodynamic properties of 2+1 dimensional, large N pure 
Yang-Mills theory living on a two-sphere of radius R. As argued in , asymptotically free 
gauge theories become weakly coupled at small volume, in this case when A = g 2 NR <C 1. 
This is because all modes with the exception of the constant mode of Aq on 5" 2 x S 1 (which 
we denote by a), are massive, with masses > 1/R. We thus have an effective IR cutoff 
on the renormalization group flow of the coupling gy M N . The massive modes may be 
integrated out directly in perturbation theory, yielding an effective action for the zero- 
mode a, or more precisely, for the constant SU(N) matrix U = e l ^ Q . The result of @ 
shows that in the large iV limit, the matrix model undergoes a Hagedorn-like transition at 
some critical temperature, which corresponds to a deconfinement transition in the gauge 
theory. We note that the large iV limit is essential for the theory on a compact manifold 
to exhibit a sharp phase transition. 
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2.1. Basic setup 

The thermal partition function of pure SU(N) Yang-Mills theory on S 2 at temperature 
T can be evaluated by the Euclidean path integral of the theory on the manifold S 2 x S 1 , 
with the radius of the circle S 1 equal to f3 = A. The action of the theory is: 

C=^j\tJ d 2 xtr(F^F^). (2.1) 

To perform this computation we will work in the gauge: 

d,A l = (2.2) 

where i = 1,2 runs over the sphere coordinates, and di are (space-time) covariant deriva- 
tives. 

This does not completely fix the gauge, as we can still make spatially independent 
gauge transformations. To completely fix the gauge we also impose 



d t A = 0. (2.3) 
Js 2 

so we choose the constant (on the sphere) mode of A to be independent of time. We 
define: 

a = ^f A , (2.4) 

cl>2 Js 2 

where uji is the volume of the 2-sphere. 

The mode a is a zero mode of the theory which is always strongly coupled and cannot 
be integrated out perturbatively. Because of this, we will compute the path integral of 
the theory in two steps. First we will integrate out all other modes which are massive to 
generate an effective action for the zero mode a. Then we will analyze the effective action 
for a. In other words, we will do the path integral in the following order: 

Z = J DAdae~ s[A ' a] = J da J dAe~ s[A > a] = J dae~ SE " [a] (2.5) 

As explained in the effective action for a can be written completely in terms of the 
unitary matrix U = e l ^ a in the form: 

S °ff =J2 C ^-m Tr (U m )Tr(U- m ) + X(3j2 C m,n,-m- n Tr(U m )Tr(U n )Tr(U- m - n )/N 
+X 2 (3 C myn ^ m . n . p Tr{U m )Tr{U n )Tr{UP)Tr{U- m - n ^)/iV 2 + . . . 

m,n,p 

(2.6) 
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Then the free energy of the theory is given by the matrix integral: 



f[dU] 



Z(J3) = e~P F = / [dU]e~ s 'ffW (2.7) 



where, as discussed in 0] , the Fadeev-Popov determinant of the gauge fixing transforms the 
integral over a to an integral over the gauge group with Haar measure [dU] . Note U is the 
holonomy of the gauge field around the temporal circle. We thus have an effective action 
for the temporal Wilson loop, whose expectation value is the standard order parameter for 
the deconfinement transition. 

In the large iV limit, we can evaluate the integral using saddle point techniques. 
Introducing the eigenvalue distribution p{6) — 5(9 — 9i)/N where Oi are the eigenvalues 
of a, i = 1..N and defining u n = J d6p(6)e inB = Tr(U n )/N the effective action takes the 
form: 

Z[(3}= J[du n }[du n }e- N2s 'eff^> u ~ n 'P^ (2.8) 

with 

&eff = ^ m Cm,—rn) u rriU'm 4" A/3 ^ ^ Cm,n, — m—n u m u n u m+n 

rn m,n ^ 

4" A (3 ^ C rn n p — m — n —pU m u ri UpU Tn j rTL j r p -\- ... 

m,n,p 

where the extra 1/m comes from the Vandermonde determinant obtained in going to the 
variables u n . 

Note that u n = is a stationary point at all temperatures. It corresponds to the 
uniform distribution of the eigenvalues of U. The stability of this saddle point depends on 
the values of the coefficients C n ^ m .... As we will show in the next subsection, the coefficient 
of |i£i| 2 is positive at small T but becomes negative as we increase the temperature. This 
signals that the u n = phase becomes unstable and the system undergoes a phase transi- 
tion. As explained in 0, near the transition temperature the u\ mode becomes massless, 
while the u n> \ modes remain massive. We will further integrate out these higher moments 
to obtain an effective action for u\ near Tjj to analyze the order of the phase transition. 
To achieve this we will need the coefficient of the quartic term in u\. The relevant terms 
in ( p.9| ) are: 

s' eff =(i - Ci,_i)K| 2 + (1/2 - c 2 ,_ 2 )M 2 + ... 

+A/3[C lili _ 2 (^u 2 + ?W) + ...] (2.10) 
+A 2 /3C 1)1 ,_ 1 ,_ 1 | Wl | 4 + ... 



5 



For the saddle point configuration, the higher modes u n are determined by minimizing the 
effective action over u n for fixed ui, so we have 

«a = -^ ii^-cl 2 ) "* + ° (A2) (2 - 11} 
This gives the effective action: 

S'effM = (1 - C^OKI 2 + A> c | Ul | 4 + ... (2.12) 



where 



6 C = -/3 2 C 1>1) _ 2 (/?) 2 /(l/2 - C 2> _ 2 (/3)) + ^Ci,i,-i,-i + 0(A) (2.13) 



As discussed in || , the sign of b evaluated at the critical temperature determines the order 
of the phase transition. If b c < the transition is first order and occurs at the Hagedorn 
temperature Th at zero coupling, but slightly below Th at small but finite coupling (they 
are the same up to order A 2 ). If b c > there are two phase transitions. The first one 
occurs at exactly Th and is second order, while the second one is a Gross- Witten type 
third order phase transition happening above Th- The third order phase transition occurs 
at the point where the eigenvalue distribution of U develops a gap. Our goal will be to 
compute the coefficients C nim ,... to the appropriate order in perturbation theory to obtain 
b c . The 6*2,-2 requires a one loop calculation, while C\^\-i and Ci^-^-i require two and 
three loop calculations respectively. 

2. 2. One loop free energy 

It is shown in [0 that the coefficient C m - m can be extracted by computing the one 
loop partition function of the theory . It turns out that: 

C m ,- m = - Z -X^l (2.14) 
m 

where x = e~$l R and zy{x) is the single particle partition function 

z v (x) = J2 n ( A ) x ~^ A ( 2 - 15 ) 

A 

In our case A 2 are the eigenvalues of the Laplacian on the unit two sphere acting on the 
vector spherical harmonics surviving the gauge fixing, while n(A) is the multiplicity of 
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each mode. More specifically = h(h + 1) , n(A h ) = 2h + 1, h = 1, 2, ... The Hagedorn 
temperature is determined by the point where the u± mode become unstable: 

z v {xc) = 1 (2-16) 
We are not able to get a closed form for zy(x), but numerically we find: 

x c ~ 1.195097 (2.17) 
so the Hagedorn temperature of the free theory on a two-sphere of radius R — 1 is: 

T H ~ 0.302675 (2.18) 

Similarly, 

C 2 ,-2 = -^p- = 0.38155 (2.19) 

2.3. Gauge fixed action 

Formally, the thermal partition function is of the form 

Z[(3] = j [da^dA'^dAijA^e- 3 ^^® (2.20) 



where Ai, A2 are the Fadeev-Popov determinants associated with (|2.2|) , (|2.3|) respectively 
and f3 = 1/T is the size of the time circle. The prime on A' signals that the constant mode 
has been removed. As shown in 0, A2 can be combined with [da] to give the integration 
measure of a unitary matrix [dU] with U = e % ^ a . On the other hand, Ai is 

detd l D l = J VcVce- tY{59iDlc) (2.21) 

where D % denotes a gauge covariant derivative 

D i c = d i -ig YM [Ai,c} (2.22) 

and c and c are complex ghosts in the adjoint representation of the gauge group. With 
our gauge choice, the quadratic terms in the Yang-Mills action ( |2.1|) take the form 

- J d 3 xtr [^Ai{bl + d 2 )A l + ^A' Q d 2 A' + cd 2 c^j (2.23) 
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where 



D X = d X-i[a,X]. 
The interaction terms in Q2.1|) are given by 

J d 3 xtr(igYMD°A i [Ai, A' ] — igyivilA 1 , A^\d\A^ — igyMdiAj[A\ A^] + 
^M^A^A^A'} - ^[A^A^A'^A^-igyMdiciA^c]). 



(2.24) 



(2.25) 



We wish to obtain the effective action by 




(2.26) 




(2.27) 



3. The Perturbative Calculation 

In this section, we summarize the computation of the two and three loop diagrams. 
As we are on S 2 , it will be convenient to expand the fields in terms of the vector and scalar 
spherical harmonics on 5" 2 . As a result, the spatial momentum integrals in the Feynman 
diagrams can be replaced by sums over the quantum numbers of the generators of SU(2). 
We will first set up the conventions that facilitate the computation. 

3.1. Spherical Harmonics Expansion on S 2 

The basic set-up for the computation was described in section 2 of 0. As in the 
3 + 1 case, it will be useful to write the action explicitly in terms of a set of spherical 
harmonic integrals. We will denote the scalar and vector spherical harmonics on S 2 by 
S a (9) and Vf(9), where a = (j a , m a ) and (3 = (jp, mp) are the SU{2) quantum numbers 
for the various modes. Our conventions for the vector spherical harmonics can be found 
in appendix A. Note for S^ ,m and V^ m the j quantum number starts at j = 0, j = 1, 
respectively. In terms of these, we have: 



a 



Ai(t,e) 



(3-1) 



c(t,9) 



J2c a (t)S«(9). 



a 
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We will denote the complex conjugate of S a by S a . On S 2 , the complex conjugation for 
scalar spherical harmonics is the same as inverting the sign of m a and multiplying by 
(— l) m ". We also denote the complex conjugate of Vf by Vf . The complex conjugation 
for vector spherical harmonics has the effect of changing the sign of mp and multiplying 
by (— l) m /3 +1 . The scalar spherical harmonics are an orthonormal basis of functions on S 2 , 



gag/3 = jotf (3.2) 

In terms of these spherical harmonics, we can define 

C apl = [ S a V p ■ VS 7 , 
Js 2 

D aPj = f v^.yPs 1 , (3.3) 
Js 2 

E a{3j = f (V X F Q ) • (fP X V 1 ) 

Js 2 

They appear as effective couplings of various interaction terms in the pure Yang Mills 
Lagrangian. These integrals can be computed explicitly using properties of SU(2). We 
will list the expressions for C, D and E in the appendix. Note that C is antisymmetric 
in a and 7, D is symmetric in a and (3, and E is antisymmetric in (3 and 7. Using these 
expressions, the quadratic part of the action for pure Yang Mills theory on the two sphere 
becomes 

S 2 = J dttr\^A*(-D 2 T +j a (j a + l))A a + ^a«j a (j a + l)a a + c^ (3.4) 

The cubic interactions are 

S 3 = 9ym J cfttr(zc"[A 7 , c p }C^ + 2ia a A^a p C a ^ (i 

- i[A a , DrA^D^ - iA a A 13 A 1 E af31 ), 
and the quartic interactions are given by 

S 4 = 9ym I dtiv{ - ha a , A^}[a\ A 6 } ( D ^ a D s ^ + , 1 - C apX CP 5X ) 

_ -A 01 APA 1 A s ^jj aiX D l3SX — D olSX D /3 ' yX ^j). 

3.2. Effective Vertices 

Since the action is quadratic in a and c, these may be integrated out explicitly to give 
additional effective vertices. As discussed in detail in M, we have two types of effective 



(3.5) 



(3.6) 
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vertices. The A type vertices arise from loops of a and c. The B type vertices are from 
open strings of a's containing two vertices linear in a and some number of quadratic 
a vertices. Both involve divergences proportional to 5(0). As for the 3+1 dimensional 
theory on S 3 , all divergent contributions proportional to 5(0) cancel. In particular the 
contributions from the a and c loops completely cancel out with the divergent parts in 
the B type effective vertices. We thus have only the non-divergent contributions from the 
latter. Computationally, we will keep only the vertices appearing in the Lagrangian above 
which contain no temporal component a and ghosts c. In addition, we will also have the 
second type of effective vertices, but can ignore any contributions proportional to 5(0) 
arising from contractions of \7 T As in these. The B type vertices are: 



The relevant diagrams contributing to the vacuum energy at one, two and three loops 
are shown in figure 2. They are the same as those in 3+1 Yang-Mills theory on S 3 . The 
B type vertices are denoted by circles. 

3.3. Propagators 

The propagators from the quadratic part of the action are 




2 ^ jx (jx + l)ja (ja + l)jr (jr + 1) ja (ja + l)jr (jr + 1) 

tr([[A a \ D r #], A 7l ][[A aa , D T A P % A 72 ]). 



(3.7) 



cUtK d (t')) 



5 af3 5(t - t')5 ad 5 ch , 



(3.8) 



a 



« b (t)a P cd (t')) = —5^5(t-t')5 ad 5 cb , 



(3.9) 



a 



(A: b (t)A" cd (t')) = 5^Af : cb (t-t')^ 
(D T AZ b (t)A p cd (t')) = -(AZ b {t)D T A? d (t>)) = 5^D T Afj cb (t - t'), 
{D T AZ b (t)D T A^ cd (t')) = 6°*6(t - t')5 ad 5 cb - 5^LlAfj cb (t - t'), 



(3.10) 



(3.11) 



(3.12) 



where the (adjoint <E> adjoint)-valued function Aj a (t) solves 



(-D 2 T+Ja (j a + l))A Ja (t) = 5 Ja (t) 



(3.13) 
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One loop: 




Two loops 






2c 



Three loops 





Fig. 2: The diagrams contributing to the free energy up to 3- loop order. In 
this figure we present a particular planar form for each diagram, but in some 
cases the same diagram may also be drawn in the plane in different ways. 

explicitly on [0, /?), 



Al d ' cb (t) = € 



ictt 



3a 



-L a t 



,L a t 



2L n V 1 - e ioL ^e- La(i 1 - e ioL ^e L ^ 



(3.14) 



where L a = (j a (j a + l)) 1 ^ 2 and a is short for a ad <g> l cb — l ad <g> a cb . The quantity that will 
appear in computation is Aj a (ti — 1 2 ) where ti and t 2 range between and (3 so that the 
argument of A Ja (t) takes value between —(3 and (3. What we will use in computation are 
the full propagators: 



A io (f) = e(f)— ( 



,(ia-L a )t 



,(ia+L a )(t-(3) 



1 



2L a 1 - e iaP e -L a fl l _ e -iaf3 e -L a f3 



)+e(-^( 



,{ia-L a ){t+(3) 



?(ia+L a )t 



2L a X l - e iafl e -L a p I _ e -iaf3 e -L a l3 

(3.15) 



) 
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2 e (iac-L a )t e (ict+L a )(t-/3) ^ e (ia-L a )(t+/3) e (ice+L a )t 

D tAj a (t) = ^(^l - ] _ e iap e -L a + l _ e -ia/3 e -L a /3) +6> ( _ ^2^~l - e i<x0 e -LaP + \ - e -ia/3 e -L„/3 ) 

(3.16) 



5.^. Two Loops 

In this section, we will compute the coefficients Ci,i,-2 from the two loop di- 
agrams. The two loop diagrams contributing to the effective potential are shown 
in figure 2. In particular, we will only need to extract the coefficients multiplying 
trUtrUtrU^ 2 + trUHrUHrU 2 in each diagram. Fortunately, these coefficients turn out 
to be non-divergent. No regularizations are required for them. The two loop diagrams can 
be computed to give 
2a: 

F 2a = Yl -^(^^ - D a ^D^)A ja (0,a ab )A jp (0,a bc ) (3 . 17) 



2b: 



'26 = 

m' s 



J2 Ml J dtAff c (t)Af ' de (t)A^' ba (t)(-E a ^E^ + 2E a ^E^) (3.18) 



2c: 

F 2c = 2_^l3g 2 YM . r -Tx (AA Ja (0, a o6 ) J D t A J> (0,a ac ) +j>(j> + l)A ia (0, « a6 )A J/3 (0, a ac )) 

(3.19) 

We have used the notation Aj a (0, a a b) to signal that the propagator participates in both 
the a and b index loop in the sense a a <8> l b — l a ® oA We isolate the coefficient of 
trUtrUtrU* 2 in each of the above and get: 
2a: 

^ (2j Q + l)(2j> + 1) e -M L ° +L f>) + e-M L «+ 2L ^ + e -^L a +L p ) 

PgyM2 h ^ 4(^(^+1)0,(^+1))^ (3 - 20) 
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2b: 



P@y_M_ J- (-A 2 (jp,j J ,j a )-2A(j p ,j 1 ,j a )A(j a ,j 1 ,jp)) 



J a >3 



^>L ot LpL 1 L a + Lp + L 1 



( -{Lf,+2LJP _j_ e -(L Q +2L /3 )/3 + e -(L^+2L Q )/3 + e -(2L a +L )f3 + ^^Lp+L^p + e -(2L 7 +L Q )/?^ 
( e -(^/3+2L 7 )/3 _ e -(L a +L 7 )/3-j ^ e -(2L Q +L 7 )/3 _ e -(L a +L )/3^ 



+— 


L a — Lp — L 7 


+ 




— L a + Lp 








-(L a +2L f3 )f3 _ e -(L /3+ L 7 ) / 3^ ) 


+ 






-(2L a +L^ 


>)0) 




— -L Q — Lp -\- Lj 




L a + Lp - 










+ 




-{L a +L )f3 _ e - 


-(2L /3 +L 1 


)/3) 



2c: 



L a + + L- 



RljjaJpJj) rh( e -(3(L a +L p ) _ e -p(L a +2L fj ) _ - 



(3.21) 



(3{2L a +L (j ) 



) 



| ^ 3p(jp + 1) jl/2( c -/3(L a +L fl ) + e -p( La +2L >3 ) + e -0(2L a +L p )^ 



(3.22) 



ja(ja + 1)' 

To obtain the above expressions, we have used several properties of the scalar and vec- 
tor spherical harmonics. These properties as well as the definitions of A(jp, j 7 , j a ), 
Ri(jp,j 1 jj a ) are listed in the appendix. We will numerically evaluate the momentum 
sums. The results are: 



Diagram 


Value 


2a 


0.000609 


2b 


-0.013068 


2c 


0.00346 



With these values we obtain: 

-/3cCi,i,- 2 (&) 2 



= -0.000816. 



(3.23) 



3.5. Three Loops 

To determine the order of the phase transition we will also need the coefficient 
Ci^-i^-i from the three loop diagrams. The three loop diagrams contributing to the 
effective potential are the same as those in 3+1 dimension. We will now list the expres- 
sions for the three loop diagrams. We will use the notation 



j^abc _ j^abc _|_ j^bca _|_ j^cab 



(3.24) 
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which is totally antisymmetric in its indicies. 

We find the following expressions for the diagrams 3a: 



r _ P9YM 

J dt A j Q ( 5 «ab)A j7 (t, a ca )A J9 (t, a ac )(A i/3 (0,a ad ) + A j/3 (0,a dc )), 

(3.25) 

3b: 



F 3b =(3g$ M {D a i x D &sx - D a&x D^ sx )'- 



3p(Jp + 1) 

{A ja (0, a ab ) (D T A Jp (0, a ad ) + D T A Jp (0, a dc )) 

J dt (D T A^(t, a ac )A js (t, a ca ) - A^(t, a ac )D T A js (t, a ca )) 

+ Aj a (0, a o6 )(A J> (0, a ad ) + A jp (0, a dc )) 

j dt {jpUp + l)A jy (t,a ac )A js (t, a ca ) - D T A^(t, a ac )D T A js (t, a co ))| 

(3.26) 



3c: 



F 3c = - 



(3g YM D a ^ x D asX D^PD s ^ 
2 j\(j\ + l)j P {j P + l) 

J dt(A jf3 (0,a ad ) + A i/3 (0, a dc )) 

{(ja(ja + l)j>0> + 1) + j' 7 (j 7 + l)js(js + l))A ja (0, a ab )A js (t, a ac )A^{t, a ca ) 

- jp(j p + l)D T A 3i (£, a ca ) {AD T A Ja (0, a a6 ) A ]s (£, a ac ) + 2A ja (0, a ab )D T A JS (£, a oc )) 
-2j 5 (j 5 + l)A ia (0, a a6 )A, 5 (0,a ac )} 

+ y dt(D T A j0 {O, a ad ) + D T A jp (0,a dc )) 

{h(h + l)A^(t,a ca )(4D T A js (t,a ac )A ja (0,a ab ) + 2A js (t, a ac )D T A ja (0,a ab )) 

- 2D T A Ja (0, a ab )D T A js (t, a ac )D T A^ (t 

i C^ca) 

-2(D T A ja (0, a ab )A js (0, a ac ) + 2A ja (0, a ab )D r A js (0, a ac ))} , 

(3.27) 
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@9yM jja^X jjpSX jja^jp jjpSp 

4 

dtA ja (t, a ab )A j/3 (t, a bc )(2A jj (t, a cd )A js (t, a da ) + A JS (t, a cd )A^ (t, a da )) 

_ @9yM jjaf3X jjjSX jjj(3p jjaSp 

4 

J dtA j a ( t ^ <Xab)Aj(,(t, «bc)(A^(t, a cd )A jg (t, a da ) - 4A js (t, a cd )A^(t, a da )), 

(3.28) 

- - - r)7<5p T~)6ipP 

F 3e =[3g YM (2D aSX D^ x - D af3x D^ x - D^D^) — 

Jp{jp + lj 

/(3 29) 
dt {D T A ja (£, a ba ) A jf3 (t, a. ad )A h (t, a cb )D T A js (t, a dc )) 

-A Ja (t, a ba )D T A jf3 (t, ot a d)A^ (t, a cb )D T A js (t, a dc ))} , 



(3g YM d^D^D^pD^p 
2~ jx(j\ + l)j P (j P + l) 
[4A Ja (0, a ab )D T A JS (0, a cd )D T A j/3 (0, a da ) 

- 2(j p (j p + 1) +j s (j 5 + l))A ja (O,a ab )A js (O,a cd )A j0 {O,a da ) 
+ J dt{A ja (t, a ab )A^ (t,a bc )A js (t,a cd )A j/3 (t,a da ) 

hih + WsUs + 1) +jf}(Jp + !)) 
+ 2D T A ja (£, a ab )D T A^ (t, a bc )D T A js (t, a cd )D T A jf3 (t, a da ) 
-Aj 1 (j 1 + l)A ja (t, a ab )A^(t, a bc )D T A j5 (t, a cd )D T A jf) {t, a da )}] (3.30) 

(3g YM D al3x D^ 6x D^PD^P 
2 j\(j\ + l)jp(j P + l) 
[4D T A Ja (0, a ah )D T A h (0, a bc )A j0 (0, a da ) 

- 2ja(j a + l)A ja (0, a ab )A i7 (0,a bc )A j(3 (0, a da ) 

- 2A Ja (0, a ab )D T A^ (0, a bc )D T A j0 (0, a da ) 

+ J dt{D T A ja (t,a ab )D T A^(t, a bc )D T A js (t, a cd )D T A j/3 (t,a da ) 

+ 3a(j a + 1)35(35 + l)(Aj a (t,a ab )A^(t, a bc )A js (t, a cd )A jf} (t, a da ) 
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+^js(js + l)(A ja (t,a ab )D T A^(t, a bc )A js (t, a cd )D T A jp (t, a da ) 
-4(j s + l) 2 D T A ja (t,a ab )D T A^(t,a bc )A js (t,a cd )A jl3 (t,a da )}] , 

3g: 

F 3g =(3g^ M E aSp E lf3p (D^ x D m - -D^ X D^ X - -D« Bx D^ x ) 

J dtdt'A j/3 (t', a da )A^(t', a cd )A js (t, a bc )A ja (t, a ab )A jp (t' - t, a ac ), 



3h: 



3i: 



3j: 



F 3 h =(3g$ M ■ , l , D^ x D^ x E^ p E^ p 
3\{j\ + 1) 

J dt 1 dt 2 D T A ja (t 1 ,a ab )A j ^(t 1 ,a bc )A jp (t 1 -t 2 ,a ca ) 

(A js (t 2 , a cd )D T A jp (t 2 , a da ) - D T A js (t 2 , a cd )A jl3 (t 2 , a da )) 

+ P9ym^-^ D af3x D^ sx E^ p E^P 

jxUx + 1) 

/**D T A J .fo,«*)A J ,(t J ,a*)A J >-t J ,a») 

{A js (t 2 ,a cd )D T A^(ti, a bc ) - D T A js (t 2 , a cd )Aj^(ti,a bc )), 

J dt 1 dt 2 dt 3 A ja (ti - t 2 , a ab )A j0 (t x - t 2 , a bc )A jy (t 3 , a cd ) 
A js (t 3 ,a da )A jp (t 1 - t 3 ,a ca )A ja {t 2 ,a ac ), 

F 3j =(5g^ M (D apx D ppx - D ppx D al3X )E^ TCT E^ f 

J dtdt'A jp (0, a ab ) A j0 (t - t', a ac )A j<y (t', a ad )A ja (t, a ca )A jT (£', a dc ), 



3k: 



F 3k =0g$ M - ( l , D^ x D^ x & pa E^ p 
J\{J\ + 1) 

J dtxdt 2 A jp (t 1 - t 2 ,a c d)A jCT (ti -t 2 ,a da ) 

{2D r A ja (0, a ab )A jp {t u a ac )D T A^(t 2 

+ A ja (0, a ab )D T A jf3 (ti, a ac )D T A^ (t 2 , a ca ) 

-ja(j a + l)A ja (0,a ab )A jf3 (t 1 ,a ac )A j „ / (t 2 ,a ca )} , 
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jjj ct f3 T j^j (3 7 p 7 a a pa r 

31 - 

J dt 1 dt 2 dt 3 A Ja (t 2 -t 3 ,a ab )A j0 (t 3 -t u a ac )A Ji (t 1 -t 2 ,a ad ) ( 3 - 36 ) 
A Jp (h , u d c) A J(T (£ 2 , «6d) A jV (£ 3 , a c& ) , 

Jja.f3\(j\apjj 1 8p^f3=y8 

F 3m , = 4(3g Y M^T~- — i iv / ■ — rrv 
JxyJx + IjJpUp + 1) 

dt(D T A ja (0, a a &) A j;3 (i, a ca ) - A ia (0, a ab )D T A jl3 (t, a ca )) 

(A j7 (t, a dc )D T A js (t, a ad ) - D T Aj^ (t, a dc )A js (t, a ad )) 

, Da5\ C \pp D ~f5pEa-f0 (3.37) 



I 



+ 2(3gYM lx(n + i)jp(j P + i) 

J dt {D T A ja (t, a ab )A jf3 (t, a hd )D T A^ (t, a da )A js (0, a ca ) 

+ 2D T A ja (t, a ab )A jfj (t, a bd )A^(t, a da )D T A js (0, a ca ) 
-Jstis + l)A ia (t, a ab )A jf3 (t, a bd )A^(t, a da )A jB (0, a ca )} 

P9YM 3 P (3 P + l)3,Ua + l)\ S h(h + l) +S hUx + l) +D 
{D T A ja (0, a cb )A h (0, a ac )D T A Jfj (0, a ad ) 

+ 2D T A ja (0, a cb )D T A^ (0, a ac ) A 3> (0, a ad ) 
-j 7 (j 7 + l)A, Q (0,a c6 )A j7 (0,a ac )A w (0,a ad )} 

^ m j^ + i)j.(j. + i) l v 3 ^TTT) +jD 

{D T A ja (0, a a6 ) Aj 7 (0, a ac )D T A Jfj (0, a ad ) 

+ 2D T Aj a (0, a o6 )D T A i7 (0, a oc ) A i/3 (0, a ad ) 
+j 7 (j 7 + l)A Ja (0,a a6 )A, 7 (0,a ac )A, /3 (0,a ad )} 

l9YM 3p(3 P + 1)3.(3, + 1) ^ h(h + 1) + 
{2D T Aj a (0, a a6 )A- Aj 7 (0, a ad ) A^ (0, a bc ) 

+(3^(3* + 1) + h(h + l))A jQ (0,« ab )A j7 (0,« ad )A J/3 (0,a 6c )} . 

(3.38) 
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Again, we will only need to extract the coefficients multiplying trUtrUtrUHrU* in each 
diagram. These coefficients turn out to be non-divergent as well. No regularizations are 
required for them. At the three loop level, the Feynman diagrams have very complicated 
expressions. However, they all have the following strucure: 

G fSim/s I fs {trU,trU^)). (3.39) 

j' s,m' s 

where Gj' Sjm ' s are the group theory factors coming from the vertices in each diagram. 
Ij' s (trU,trU') come from the propagators. We will expand Ij> s (trU,trU^) in powers of 
trU, trU^ to extract the coefficients of trUtrUtrU'trll* . Fortunately, the relevant terms 
m are very similar to those in 3+1 dimensions, which have been computed 

0. To apply them, we only need to replace the masses of the propagators (j a + l)§+i — » 
(ja(ja + l))2+i- With these we can evaluate the angular momentum sums numerically. 
The results are 



Diagram 


Value 


Diagram 


Value 


3a 


-0.000182 


3b 


0.000008 


3c 


0.000432 


3d 


-0.000187 


3e 


-0.0001396 


3/ 


0.0001399 


3<7 


0.0021095 


3h 


0.00685 


3i 


-0.005715 


3j 


0.00197 


3k 


0.000779 


3/ 


-0.00195 


3m 


0.000752 


3n 


0.001682 



The total three loop contribution is 



C h i _i _i = 0.001237 (3.40) 

We thus obtain the coefficient b c = 0.00503 > 0. Thus, we conclude that the deconfine- 
ment transition is second order, followed at a slightly higher temperature by a continuous 
transition in which the eigenvalue distribution develops a gap. 

4. High Temperature Limit 

In this section, we consider the behavior of the theory at general values of the spatial 
radius in the limit where the temperature is very large (compared with either the inverse 
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spatial radius or the gauge coupling). In this limit, the Kaluza-Klein modes on the thermal 
circle become very massive, and the theory is well described by an effective two-dimensional 
theory on S 2 . This theory contains a two-dimensional gauge field together with an adjoint 
scalar coming from the zero mode of Aq on the thermal circle. As in the more familiar 3+1 
dimensional case, this scalar receives a Debye mass at one loop, corresponding to screening 
of electric charge. In our case, the mass was calculated in [|J to be 

m 2 ~ AT In 

where the logarithm arises from a resummation of infrared divergent diagram sHI. For T > 
A, this mass is much larger than the scale M ~ V^2 = V AT associated with the two- 
dimensional Yang-Mills theory, so the model should be effectively described by pure two- 
dimensional Yang-Mills theory on S 2 (as long as the sphere radius is larger than m _1 ). 

_ i 

This theory has a third-order phase transition at R ~ A 2 2 m _1 [Q, so we conclude that 
our phase diagram has an additional phase transition line coming from large temperature 
along the curve 

TR ~ 1/(XR) . 

An important question is whether this transition remains sharp for large but finite values 
of the temperature or whether it becomes smoothed out, either for any non-infinite value 
of the temperature (i.e. any finite Aq mass) or below some particular temperature. For 
high temperatures, the question should correspond to asking about the fate of the phase 
transition in the two-dimensional theory when the mass of an adjoint scalar is reduced from 
infinity. Unfortunately, even this question seems difficult to approach since the theory is 
no longer solvable with the adjoint scalar. 

The persistence of a sharp phase transition would be guaranteed if there were some 
order parameter associated with the transition. In the pure two-dimensional Yang-Mills 
theory on S 2 , there is in some sense an order parameter for the phase transition, namely 
the eigenvalue distribution for a maximal area Wilson loop which divides the sphere into 

2 Strictly speaking, there are no infrared divergences since we are working on a spatial sphere. 
Nevertheless, this infinite volume result should be valid as long as the geometrical infrared cutoff 
scale 1/R is smaller than the dynamical infrared cutoff scale m. In this case, we are still required 
to resum a large number of (finite) diagrams to get the leading contribution to the mass. For 
smaller radii, the geometrical infrared cutoff dominates and the effective scalar mass will be given 
by the one-loop contribution. 




19 



two equal areas.! This eigenvalue distribution is gapped for radii less than the critical 
radius, but ungapped above it @]. Away from infinite temperature (where the full theory 
is on S 2 x S 1 ), it will presumably no longer be true that the expectation value of a spatial 
Wilson loop will depend only on the enclosed area, so the extension of the order parameter 
to finite temperatures is ambiguous. We could for example choose to focus on the Wilson 
loop around the equator of the S 2 , and it will certainly be true that the full phase diagram 
will divide into regions for which the eigenvalue distribution for this Wilson loop is gapped 
or ungapped, but it isn't clear that the boundary of this region should correspond to some 
non-smooth behavior of the free energy away from the infinite temperature limit. 

If the high-temperature phase transition does remain sharp, an intriguing possibility 
is that it connects on to the third order phase transition line that originates at the critical 
temperature in the small volume limit (as in figure 6b). Indeed, in the limits where we 
have analytic control, both of these phase transitions are third-order transitions associated 
with gapping for the eigenvalue distribution of Wilson lines. To investigate whether the 
two transitions might be the same, one approach would be to study the behavior of the 
eigenvalue distribution for the Polyakov loop in the vicinity of the high-temperature transi- 
tion. This distribution should be close to a delta function in this high-temperature regime, 
but could still be either gapped or ungapped depending on how the eigenvalue distribution 
falls off away from the peak. A transition between these two possibilities may show up as 
a change in behavior for the eigenvalue distribution of the massive adjoint scalar in the 
effective two-dimensional theory, for example from a strictly localized distribution to one 
with an exponential tail. Unfortunately, we have not been able to determine whether such 
a transition occurs, so we leave it as a question for future work to determine whether the 
two third-order transitions are connected. 

5. Possible Phase Diagrams 

We have seen that for small sphere volumes, our gauge theory undergoes a second order 
deconfinement transition followed by a third order gapping transition as the temperature 
is increased. We would now like to understand in general the simplest possibilities for what 
happens to this behavior as the spatial volume is increased. 

3 The precise shape of the loop is unimportant since the theory is invariant under area- 
preserving diffeomorphisms. 



20 



First, we expect of course that the deconfinement transition extends all the way to 
large volume, where the transition temperature should be of order A. The simplest log- 
ical possibility is that the qualitative behavior we found is unchanged as we go to large 
volume. This would mean a second order deconfinement transition at large volume. By 
the conjecture of Svetitsky and Yaffe ||, the critical behavior should then be the same 
as a two-dimensional spin model invariant under the same global symmetry, in this case 
Zn^oo ~ U(l). This corresponds to an XY model, for which the transition should be of 
Kosterlitz-Thouless type. 

However, if the expected extrapolation to large iV of the lattice results mentioned in 
the introduction holds, the infinite volume transition should be first order for large N, so 
the Svetitsky- Yaffe predictions would not apply. In this case, which we will assume for the 
remainder of the section, there must be a critical sphere radius at which the deconfinement 
transition changes from second order to first order. 

We will now argue that there are two different types of behavior possible at the critical 
radius. To understand this, consider the effective potential for the eigenvalue distribution 
evaluated at the deconfinement transition temperature for various values of the radius. 
Where the transition is second order, the potential at the transition temperature has a 
global minimum for the uniform eigenvalue distribution, with a vanishing second derivative 
along some direction. The effective potential develops a negative second derivative along 
this direction as we go above the transition temperature, and the global minimum smoothly 
moves away from the uniform eigenvalue distribution. There are two qualitatively different 
effects that can give rise to a change to first order behavior: 

First, the effective potential evaluated at the transition temperature might develop a 
second global minimum at a point away from the uniform distribution for some value of 
the radius. In this case, if we move further along the line where the uniform distribution is 
marginally stable (dotted line in figure 3), this new minimum will (generically) become the 
global minimum, so we must have a first order transition to this new minimum occurring 
at some lower temperature. Thus, below this critical radius, the deconfinement transition 
is second order and follows the boundary along which a local instability develops around 
the uniform distribution. Above the critical radius, the transition is first order and follows 
the boundary along which we have two global minima. Generically this line of two minima 
will not terminate at the critical radius but will continue to smaller radii above the decon- 
finement temperature. Even at these smaller radii, it represents a phase transition, since 
below the line the minimum near the uniform distribution should be the global minimum, 
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while above the line, the other minimum will be the global minimum. Thus, we have a 
triple point at the critical temperature and critical radius separating three distinct phases. 
If the second minimum is at the boundary of configuration space, the higher temperature 
transition will correspond to a gapping transition, but in this case a first order one, since 
the eigenvalue distribution jumps discontinuously. The phase diagram in the vicinity of 
the triple point is sketched in figure 3. 




Fig. 3: Phase diagram in the vicinity of the deconfinement transition when 
we have a triple point at the critical radius, with sketches of the effective po- 
tential in each region. Deconfinement transition switches from second order 
(dashed line) to first order (solid line). Dotted line is not a phase transi- 
tion but represents boundary in deconfined phase of region for which local 
minimum exists at the origin. 

We now explain the other possible behavior near the critical radius (depicted in figure 
4). If we follow the curve along which a single marginally stable direction exists (with all 
other directions stable) , it may happen that this direction becomes marginally unstable at 
some point (e.g. if the fourth derivative of the potential in the marginal direction switches 
from positive to negative). Such a point is known as a tricritical point. In this case, if we 
continue along the line where we have a single marginal direction, the global minimum of 
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the potential will shift gradually away from the uniform distribution. The deconfinement 
transition then no longer coincides with the curve along which the uniform distribution 
becomes locally unstable, but occurs at some lower temperature where the minimum at 
the origin (which exists everywhere below the marginal stability line) takes the same value 
as the new nearby minimum. The deconfinement transition corresponds to a jump from 
the minimum at the origin (i.e. the uniform distribution) to the nearby minimum. It is 
therefore first order but with a latent heat which vanishes as we approach the tricritical 
point, where the two minima in the effective potential merge. Unlike the other scenario, 
there is no other phase boundary emerging from the point at which the deconfinement 
transition switches behavior, since the line along which we have two equivalent minima 
simply ends at the tricritical point. In the present case, the first order deconfinement 
transition either to the left or right of the tricritical point is certainly to an ungapped 
phase. 
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Fig. 4: Phase diagram near the critical radius in the case when we have 
a tricritical point, with sketch of the effective potential in each region. De- 
confinement transition switches from second order (dashed line) to first order 
(solid line). Dotted line is not a phase transition but represents boundary in 
deconfined phase of region for which local minimum exists at the origin. 
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Both of these scenarios are realized in a simple toy model for a complex order param- 
eter w with effective potential 



S e ff = a\w\ 2 + b\w\ 4 + c\w\ 6 

and a boundary \w\ = 1 for the configuration space. The phase diagram for this toy model 
as a function of the parameters a, b is shown in figure 5 for the two cases c < and c > 0. 
In the first case, the switch from second order to first order behavior for the transition 
corresponds to a triple point, while in the second case, it corresponds to a tricritical point. 
Note that in the case c > 0, only the details of the potential near w = are important, 
so any higher order terms can be ignored (so long as they do not give rise to a lower 
minimum). Thus, the behavior of the toy model in the vicinity of the tricritical point 
should precisely coincide with the behavior of the Yang-Mills theory if we have a tricritical 
point at the critical radius, since the effective action for u\ will be of this form, with higher 
order terms that we can ignore. 



n 



III 



c<0 



c>0 



Fig. 5: Phase diagram for toy model effective potential for c < and 
c > exhibiting the triple point and tricritical point behaviors. Phases I, 
II, and III correspond to having the global minimum at the origin, in the 
bulk of the configuration space away from the origin, and at the boundary of 
configuration space respectively. Solid and dashed lines represent first order 
and second order phase transitions respectively. 
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We can now comment on the possible forms for the full phase diagram, assuming that 
the large- volume transition is first order, so that there exists a change of behavior at some 
critical radius. If the effective potential is such that we have a triple point, the simplest 
possibility would be that the additional phase boundary coming from the tricritical point 
corresponds to a gapping transition, so that this phase boundary would connect with the 
third order gapping transition emerging from the zero-volume critical temperature. This 
requires there to be some radius at which the gapping transition switches from third order 
to first order. In this scenario, the full phase diagram would appear as in figure 6a. The 
fate of the high-temperature Douglas-Kazakov transition is unclear, though as we have 
discussed, it is possible that this phase boundary simply ends at some high temperature. 

In the case where the effective potential gives rise to a tricritical point, the gapping 
phase boundary could either extend up to infinite temperature and be absent in the large 
volume limit (e.g. in the scenario where it connects with the Douglas-Kazakov transition 
at high temperatures), extend to large volume such that the first order deconfinement 
transition there would be followed by a gapping transition at some higher temperature, 
or end somewhere on the deconfinement phase boundary to the right of the tricritical 
point (as in the toy model for c > 0). These possibilities are shown in figure 6 b,c, and d 
respectively. 

Further lattice studies should help distinguish between the possibilities in figure 6. In 
particular, while the distinction between gapped and ungapped eigenvalue distributions 
strictly exists only in the large N limit, recent studies at relatively large but finite values 
of N have provided clear suggestions of gapping transitions for eigenvalue distributions 
of spatial Wilson loops (see, for example |1(J). Thus, it should be possible to determine 
whether the deconfinement transition at large volume is to a gapped or ungapped eigen- 
value distribution, and in the latter case, whether there is an additional gapping transition 
at higher temperature (as in figure 6c). On the other hand, distinguishing between pos- 
sibilities a) and d) may be difficult, since they differ only at intermediate values of the 
radius / coupling. 

6. Conclusions 

The main result of this paper is that pure large N two-dimensional Yang-Mills theory 
has a second order deconfinement transition at small spatial volume, with a third-order 
gapping transition at some higher temperature. This is is a qualitatively different behavior 
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Fig. 6: Simplest possible phase diagrams for large N pure Yang-Mills theory 
on S 2 as a function of sphere radius R and temperature T, assuming a first 
order deconfinement transition at large volume. Solid, dashed, and dotted 
lines correspond to first, second, and third order transitions respectively. 

from pure Yang-Mills theory in 3+1 dimensions, and provides the first example of a gauge 
theory with a single-trace Lagrangian in more than two space-time dimensions for which 
the deconfinement transition is second order at small volume, and which therefore displays 
three distinct phases. 

If the same behavior exists for some large N theory with a controllable gravity dual, 
it would be fascinating to understand what the new ungapped phase corresponds to on 
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the gravity side. Generally, the deconfined phase of a large N gauge theory corresponds 
to a black hole geometry.0 In J\f = 4 SYM theory at strong coupling, the deconfinement 
transition corresponds on the gravity side to a first-order transition between the original 
AdS 5 x S 5 spacetime with a thermal gas of supergravity particles to a large black hole 



spacetime ||11|| . On the other hand, in a theory with an intermediate ungapped phase on 
the field theory side, there should be a stable intermediate type of black-hole phase on the 
gravity side smoothly connected to both the no black hole phase and the big black hole 
phase. 
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Appendix A. Spherical Harmonics 



In this section we set up our conventions for the vector spherical harmonics , . 
The definition of a vector spherical harmonics is: 

Vjim = V JiM^i = Y i™e g (l mlq\JM) (A.l) 

q m,q 

where the e q are in the spherical tensor basis: 

+ ie y 



2V2 



2V2 

e = e z . 



(A.2) 



4 The argument Qis that the non-zero expectation value for the Polyakov loop implies that a 
string worldsheet whose boundary wraps the thermal circle in the associated Euclidean spacetime 
can have finite area, and therefore the thermal circle must be contractible. In the Lorentzian 
picture, this is associated with the existence of a horizon. 
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The raising and lowering of the vector index q is given by: 



Vjim, 9 = (-i) 9 ^7m (A.3) 

and 

Vj lM = {-l) M+J ~ l+1 Vji-M. (A.4) 

The vector spherical harmonics can be used to expand any well behaved vector fields in 
R 3 and they can be categorized into the following orthonormal basis [ 13| : 



PjM = (2J l l)1/2 [~( J + !) 1/2 K/J+1M + J 1/2 Vjj- 1M } = fYjM 



(2J + 

Bjm = , OT , lU/2 [ ^ ' V JJ+1M + (J + 1) Vjj_ 1M \ = VYj M (A.5) 

(2J + 1) L / 2 r r 

r - 

Cjm = — r x — — /9 Bjm = —iVjjM- 

(J(J + I)) 1 / 2 

The first of the above does not live in the tangent space of the two sphere, while the 
second do not contribute to the effective action by gauge fixing. We have the following 
useful expression for Vjjm following from the last of the above identities: 

V JJM,q = + \}y/2 YjM +i ( A - 6 ) 

1 // 7 7 i ~\jt i i\\l/2 _ ( _1 ) J_M //o T , Q wo 7 , i\fn 7^^l/2 ( J J 1 

M -M - 1 1 



— 1 (—~\\ J ~ M ( 
L + = -^((J-M)(J + M + l))Va = _LJ1 ((2J + 2)(2J + 1)(2 J)) 1 / 2 ^ 

1 f_lU+M / 7 7 1 

^- = ^T72(U + M)(J-M + l)) 1 /2 = ^ ((2J + 2)(2J + l)(2J)) 1 /^_^ ^_ j J 

L = M = ( ~ 1) 2 J " M ((2 J + 2)(2 J + 1)(2J)) 1 / 2 ( ^ _ J M 

(A.7) 



Appendix B. Effective Vertices 

As we have seen in section 3, when expended in terms of scalar and vector spherical 
harmonics, the Lagrangian for pure Yang-Mills theory on S 2 x S 1 contains effective vertices 
that are integrals of products of spherical harmonics. Here we explicitly compute effective 
vertices. We will write the results with 3-j symbols: 

(i 1 i 2 i 3 )=(-l) Jl - J2 - m n^3 + l)- 1/2 (jirn 1 j2rn 2 \j 3 -m 3 ). (B.l) 
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A useful formula for sum of products of three 3-j symbols: 

h h h \ ( h h h 



m 2 /i 3 



h h h \ _ ( h 32 33 \ (ji h h 
Mi -A*2 rn 3 J \rrit m 2 m 3 J \ h h h 

SSS (given by Gaunt's formula for associated Legendre polynomials) fl 



x 



(B.2) 



v v v _, (2Zi + l)(2f 2 + l)(2Z 3 + l) u/2 ^ 1 Z 2 h\(h h h\ 



S 2 

h h h 



zI hhi 3 



mi 7712 m 3 



(B.3) 



VSV, The D vertex 



^ aiQ2Q3 = / V hhmi V l2l2m2 Y hm3 = 
Js 2 

(2h + 1)(2Z 2 + 1)(2Z 3 + lh i/a jagg + jj ~ flgi + !) ~ jaCjg + 1) 

4tt J 2(Z 1 (Z 1 + 1)Z 2 (Z 2 + 1))V2 ^ 

/i k h\(h h l A=R\/\ h ^ h )( 1 ^ l * l * 
(J (J J V mi m 2 m 3 / ymi m 2 m 3 

Note Zi + / 2 + Z 3 has to be even for nonzero amplitudes. 

• SVS, The C vertex 



C a3a2ai = [ (VY hmi )V hhm2 Y hm] 
Js 2 



1 (2Zi + 1)(2Z 2 + 1)(2Z 3 + l) )1/2( (J + 1)( J - 2Z 3 )(J - 2Z 2 )(J - 2Zj + 1) 1/2 
2r l 4tt J 1 Z 2 (Z 2 + 1) ' 

Zi-1 Z 2 Z 3 \ / Zi z 2 z 3 \ =A f h h h\ 
J I mi m 2 m 3 y 123 V mi m 2 m 3 y 

(B.5) 

where J = l\ + Z 2 + Z 3 and it has to be odd for non-zero amplitudes. 
• (curlV-r)(VxV-r), The E vertex 

E « 3 aia 2 = /" ( V X Vj 3 J 3 M 3 • r)(K/ lJlMl X V> a j aMa • f) 

_ J 3 (J 3 + 1) 1/2 / Ji J 2 J 3 \ 

- _( Ji(Ji + l) ) AjlJ2j3 l,Mi M 2 M 3 J 

T ( J\ J2 Jz 

Note since J = Ji + J 2 + J 3 has to be odd in order for Aj x j 2 j 3 to be nonzero, 

fl ft fl\ wm P^k up a negative sign upon interchanging J±, J 2 . This suggests 
IVli 1VI2 M\\ I 

Aj 1 j 2 j 3 needs to be symmetric in Ji, J 2 , which can be checked to be true in our expression. 
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Appendix C. Summation Formulas 

The following identities are useful in computing the two loop diagrams. 

J2 Da ^ D ^ = Mja,jp,j^) (CI) 



47T 



where 



„ , , , , (2h + l)(2/ 2 + 1)(2Z 3 + 1) , (h(h + 1) ~ h(h + 1) ~ h{k + I)) 2 
* l(/l ' /2 ' /3) - ( ^ } 4/ l( / 1 + l)/ 2 (/ 2 + l) 

J 

Y^ E ^ E ^i = +A 2 (j p ,j 1 , Ja )a(J a ,J p ,J J ) (C.5) 

m's 

Y^E^E*** = (-l)A(j P , J 7 , ja)A(j fll J 7 , J»Or(J a , ^, J 7 ) (C.6) 
m's 

where /2, h) is defined as above, and cr( J a , Jq, J 7 ), is 1 if its arguments satisfies the 
triangle inequality, and is zero otherwise. The presence of cr( J a , Jg, J 7 ) is just a reminder 
that we need to impose the triangle inequality on J a , Jp, J 7 at each vertex which is obvious 
from the left hand side but is somewhat obscured by the summation. 

The following identities are useful in computing the three loop diagrams. Again, a 
hat over the indicies means summing over their cyclic permutations. 

££>* 3A = -(-l) J «flJ/ 2 (J«, J a ,0)(2J a + ly/Hj^Sm^o (c.7) 

m 

J^D^D^ = (-ly^ —L— Sj^jMa,^) (c.8) 

m's V 7 

E ^ a7T ^~ = (-1) J " +1 (2 1} Jj a , J,fli(a, 7, r) (C.9) 

m's a 

E^ 7 ^ = (-1) J " +1 1) ^ a ,J^( 7 ea)A(a6 7 ) (CIO) 

m's V a 
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Y^E^E^ = (~1) J " (2J * + Sj^A^ea^ae) (C .H) 

m's a 

E D7 ™^ = ° (C.12) 

m's 

E^ 7T ^ = ° (C.13) 

m's 

jjaaX jj-yeX jj-yer jjf3(3r 

m' s,J\,J T 

= (_i)^+^^ Je ( 2 J Q + i)i/2(2 J/3 + l^R^ia, a, 0)< 2 (/3, 0, 0)fli( 7 , 7, 0) 
= 5j 7 ,j eT ^(2J Q + 1)(2J 7 + l)(2J /3 + 1) 

(C.14) 

J- D aaX W eX D ^r D P^ = (_!) ^^^ V «, 0)< 2 ( 7 , 7 , 0)^(0, 7 , r) 

m's,J A ^ + iJ 

(C.15) 

Y^D^D^D^D^ = 6j^j. 1 + fli(a, 7 , A)fli(^ 7 ,r) ( C .16) 

m's V 7 

YD^D^D^D 7 ^ = 8 Jx ,j T — ^—rR 1 (a, 1 ,X)Ri(e,P,X) (C.17) 

m's 

^ D af3X D jeX D Ter D aW = _{£ jt s R l/2 ^ ^ X)R[ /2 (-y, e? A)#} /2 (/3, 6, r)i?J /2 (a, 7 , t) 
m's 

I «Xy i/q/ 1/7- I 

\ J/3 J e Jx / 

(C.18) 

E D^D^E^'E^ = -(-l) Ja+Jp Sj a ,j p f Jp ± f'l R{ /2 (a, a, 0)< 2 (p, p, 0) 

m's,J A [AJa + ^ 

A — A — 

TOOL UTOL 

(C.19) 
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Y^C^C^D^D^ = d Ja ,j ^ - V^i(«,^,r) (C.20) 

m's T 

^^a^ 7 a^p^ = _^ >Jfl 3_^_ fll(a)7)A) (a21) 

m's V 7 

J-,,Jz 2T + 1) «P7 aep r/3 7 re/3 ^.ZZj 
m's V 7 ' 

m's L e Z 3 T J 

(C.23) 

V E paT E^E^EP a ° = -lEi'.A^ A-^A^A — { J T P J T a J T T 1 fc 24) 



£ DorrppfhacpfiTCTa* = #1/2^ % p)R \/^^ 7> a)A Tf3p A aaT I ^ £ ^ } 

m's I ^ T CT 

(C.25) 



Y^E a ^C /r " 7 D 0T<T D aT " = -{-l)^i'°A-^A alp R\ /2 {p, r,a)i?; /2 (a,r,p) | 



j 

- Jaj 

(C.26) 



expressions that are zero : 
in 3h: 

Y^D^D^E^E^ = 5 Jx ,j * + < 2 (a, 7 , A)i?} /2 (/?, 6, X)A^A^ = 

m's 

(C.27) 

in 3m, 3g: 

£ jjafiX^jyyep^ = ^ _J_ fl l/2 (a> p? A )^ /2 ( 7 , 6, p)A paA 2^ = (C.28) 



We can use the following to simplify the above expressions 

8tt 



EiMM,W= (2 ' 1 + 1 ' (2 ' 2 + 1) (C.29) 



*3 



< 2 (/ 1 ,/ 1 ,0) = (-l)^ l (2 ^+J 2 )1/2 (C.30) 
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